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INTERVAL-CENSORED  TYPE  H  PLAN 


Surekha  Mudivarthy,  M.  Bharkara  Rao,  North  Dakota  State  University 
Rupa  Mitra,  Moorhead  State  University 

MLB.  Rao,  North  Dakota  State  University,  Department  of  Statistics,  P.O.  Box  5575,  Fargo,  ND  58105 


Key  Words:  Censored  data,  efficiency,  interval- 
censored  data.  Type  II  plan 

1.  Introduction:  The  Type  II  plan  is  used  on  many 
occasions  in  order  to  estimate  the  lifetime 
distribution  of  a  product  under  investigation.  Begin 
the  process  by  choosing  two  positive  integers  r  <  n. 

Select  a  sample  of  n  units  of  the  product,  set  them  to 
work,  and  observe  the  units  continuously  until  r 
units  fail.  The  objective  is  to  estimate  the  lifetime 
distribution  of  the  product  using  the  data  on  r  failure 
times. 

We  want  to  offer  a  modification  of  this  plan  in 
response  to  a  past  consultation  problem.  This 
problem  arose  from  two  diverse  fields:  one  from 
engineering  and  the  other  from  ornithology.  In 
order  to  expedite  observation  of  the  r  units  failure 
times,  only  periodic  inspections  were  made.  In  such 
a  plan,  the  exact  failure  times  of  the  r  units  will  be 
unknown,  but  we  will  know  how  many  units  failed 
between  each  of  the  consecutive  inspection  times. 

Formally,  the  inspection  planxan  be  described 
as  follows.  Choose  and  fix  a  number  /q  >  0  .  Select 
a  sample  of  n  units  and  set  them  to  work.  Inspect 
the  units  at  times  (q  ,2/q  , . . .  until  r  units  fail.  Let 
M  denote  the  number  of  inspections  needed.  Let 
=  Number  of  units  failed  during  the  i^  inspection 

interval  ((/ - 

/=1,2,3,... .  The  data  consist  of 

...,Xj^.:r  - 

These  random  variables  satisfy  the  following 
conditions: 

Xi  +  X2+  .  +X^_i  <^ - 1, 
and 

X]  +  Xi+.-.+Xf^^  >r. 

We  call  this  plan  Interval-Censored  Type  II  plan. 

The  basic  questions  we  were  asked  to  address  were: 

1.  Evaluate  the  loss  of  information  in  some 
meaningful  way  if  one  adopts  the  Interval-Censored 
Type  II  plan  over  the  traditional  Type  II  plan; 

2.  Provide  some  guidelines  as  to  the  choice  of  Iq  . 


In  this  paper,  we  attempt  to  answer  these 
questions  to  the  best  of  our  ability.  This  is  ongoing 
work  and  we  want  to  report  what  we  have  achieved 
so  far.  Let  Tbe  the  underlying  lifetime  variable 
associated  with  the  product.  We  assume  that  T  has 
an  exponential  distribution  with  unknown  parameter 
^  >  0 .  The  probability  density  Ametion  of  T  is  given 
by 

f0{x)  =  ee~^ ,  x>o,  e>o. 

We  obtain  the  likelihood  estimator  of  6  under  the 
Interval-Censored  Type  II  plan  and  compare  its 
variance  with  the  variance  of  the  unbiased  estimator 
of  0  built  on  the  maximum  likelihood  estimator  of 
0  under  the  Type  II  plan.  This  comparison  is  done 
using  extensive  simulation  studies.  Some  relevant 
distribution  theory  and  asymptotics  are  presented 
here. 

2.  Some  Distribution  Theory:  Let 
be  ltd  copies  of  T  and  ^1)  <  ^2)  ^(n)  ‘*'6 

corresponding  order  statistics.  Under  the  Type  II 
plan,  the  data  consist  of  7^i) ,  7^2)  >  •  ■ » 
maximum  likelihood  estimator  of  0  is  given  by 

Under  6^  2rOI0^^  has  a  chi-squared  distribution 

with  2r  degree  of  freedom.  The  estimator  0,.  ^  is 
biased,  but  the  bias  is  correctable.  More  precisely, 

=  (r/(r  -  1))^  for  all  e>0. 

The  unbiased  estimator  ((r  -  1)/^)^/-^  has 

variance  Jir  -  2)  provided  r  >2  .  (Epstein 
and  Sobel  (1953).) 

Let  us  focus  on  the  Interval-Censored  Type 
II  plan.  We  need  to  determine  the  joint  distribution 
of  the  data  M,  Xi,  X2,  X^^,  We  proceed  in 

two  stages.  First,  we  obtain  the  distribution  of  A/ 
and  then  the  conditional  distribution  of 

x,,.Y2, 


-  a 


Distribution  of  M:  Pr^(A/  =  1)  =  Pr^(Xi>  r) 
=  Ftff  (at  least  r  failures  in  the  interval  (0,  (q  ]) 


P\ 

=  J 


-z'-\\-zy-''dz 


0  B(r,n-r  +  \) 

=  Ip^(.r,n-r  +  \), 

where  /?j  =  and  is  the 

Incomplete  Beta  Function.  (Abramowitz  and  Stegun 
(1965).)  For  m>2, 

Pr^CA^  =  nt) 

=  Pr^CA'i  +  X2+...+X^_^<r  -  I 
and  X^  +  X2  /■) 

r-1  n-s 

=  Z  Z  Pr^  utiits  fail  in  the  interval 

5=0  f=r-5 

(0,  (/w  -  l)lo  ]  and  t  units  fail  in  the  interval 
r-1  n-s 

=  z  z 


(■ 


al _ 

-.?-/)!  V  / 


n\ 


r-1 

5?o 


r-1 

j?0  5!(w-5)! 


M_^l_^-^(m-l)ro_y 

(£-5)! 
n- 
«! 


s  '\n-s 

Pm-\{^-Pm-l) 


n-s-t 


{r-s,n-r  +  \). 


where  /?,•  =^l-e  /=1,2,... . 


Conditional  distribution  of 
Under  t?>0, 

?X0iX^=x^,X2=X2,...,X„=x„\M  =  m) 


n\ 

Xi!x2!...x„!(w-Zr=iJC,)! 

^l  e  1 

^e'^o  _  ^g-(.m-l)to0  _^-mto0y’- 

forall  0<Xi,  X2,  ...,  x„<n,  2  Xj<r-l,  and 

;  =  1 
m 

X  Xi>r.  For  the  sake  of  simplicity,  let 
/=!  ” 

m 

Xm+l  =«-Z^/  • 

1  =  1 

We  now  derive  the  maximum  likelihood 
estimate  of  0  based  on  the  data:  number  of 
inspections  made  and  the  number  of  units  failed  in 
each  inspection  interval,  i.e., 

M  =  m,A'i  =x,,A'2  =X2,...,X„  =x„. 

The  likelihood  L  of  the  data  is: 

L  =  Constant(l-e'^»)'''  (e~^° 


-(m-l)fo0  _p-mtQ0 


m+l 


JC,+X2+...+JC„ 


=  Constant^l  -  j 

+2.1C3  +3x4  +...+mx^t 

The  log  likelihood  is  given  by: 

In  L  =  Constant  +  ^  Z  j  j 

-(6lfo)[  z«/+i]  • 

The  derivative  of  the  log  likelihood  is  set  equal  to 
zero  in  order  to  obtain  the  maximum  likelihood 
estimate.  The  following  equations  achieve  the 
objective. 


dd 


(lnl)=  Z^,- 
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-^1  Z«,>i  I  =  0 

</=l 

m 

z«,+i 

_i=i _ 


m 

'Lxi 
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S^+i 

-OtQ  _  »=1 
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^=-— Inl 
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m  m 

ZXj  +  I/X,+, 

1=1  1=1 

r  m 

TiXi+i 
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—In 


m  m 

I.Xi  +  S/x,-+i 

'^1=1  1=1  ^ 

Z'  m  \ 

T.Xi 
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m 

,  S«,+1 , 
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The  maximum  likelihood  estimator  of  6 ,  upon 
replacing  the  data  by  the  corresponding  random 
variables,  is  given  by 

/'  M  \ 


1 

^  =  — In 


1  + 


/=! 


M 

,  Z'Xm 

/=! 


The  next  objective  is  to  obtain  the  asymptotic 
variance  of  d ,  Rewrite  the  derivative  of  the  log 
likelihood  as 


^(lnI)  =  fo  ZXA[ - 

36  V/=i  Ai-e  “o  . 


m  ^ 

-^1  Z^^j+i 
^/=i  y 


from  which  we  have 

A(i„i)  =  -,5|  ^x. 


ae 


fK4 

Z- 

M=1 


The  asymptotic  variance  of  0  is  given  by 

Asyvar^(^)  =  —  ^ 


(  ^  ^ 
Ee  -^(InZ) 


36^ 


~~  jn 

tlEel 

^i=\ 


The  formula  for  the  asymptotic  variance  simplifies  to 
evaluating  successfully  A',  j .  We  will 

evaluate  the  expectation  using  the  conditional 
expectation  argument.  Note  that 

(M^  (  f  M  ,  \\ 

Ee\-LXA=I^E 


V|=l 


V/=1 


JJ 


m 


=  ZE\  =  Pr^(M  =  /n), 

m>\  v/=l  / 

The  critical  step  is  the  evaluation  of  the  conditional 
expectation; 


=  Z(xi  +X2+...+X^)- 


n\ 


x,!x2!...x^!x^+,! 

^e-^o  _  g-2toOy\  ^g-('»-i)toO  _  g~mtoOy- 
where  the  summation  is  taken  over  all  integers 

>  <^2  >  •  *  *  >  y  ^m+\  ~ 
Xi,X2,...,x„+i<r-l 
Xi,X2,...,x„>r, 
and 

Writing  ,-^2  >  ^  ^ »  we  can 

rewrite  the  conditional  expectation  as 

^[z^/|A/  =  /nj 

r-l  n-s  m| 

=  Z  Us+t) 


,=o  t=r-s  s\t\(n-s-t)\ 

s\ 


X, +X2 


m-l  ■ 


L-(«-l)'o«  _  L-”<0»V 


r-l  n-s 

=  Z  Z(^+0 


n\ 


s=0t=r-s  s\t\(n-S-t)\ 


-  e y  [e “'"^0^  J 


n-s-t 


n\ 


r-1  ns 

=  Z  Z{s+t)- 

s=0t=r-s  Slj!(/7-S-t)! 

^1  _  g-(n>-i)toffy  ^g-(m-i)/o0  _  g-'^fo^y 
^g-mtoOy-^-‘ 

=  sV) - - - (i-e-<'”-')'o<^ 

“o,=7-5  5!/!(«-^-/)!^ 


n-s-t 


i=0/=r-5  5!/!(«-i-/)! 


'■-1  (w-1)! 


ns-t 


V  f.-r\(x  — 


5=o(^-1)!(«-0! 

■^“1  /  /_  i\» 


Ip^(r-s,n-r  +  l) 

+(n  -  5)(l  -  )  Y  - ~  - • 

(i-e-^oy-‘(,-«o)"- 


\n-s-t 


estimate  of  0  computed,  and  the  empirical  mean 
square  error  evaluated.  The  objective  is  to  compare 
the  mean  square  with  the  variance  of  the  unbiased 
estimator  of  6  based  on  the  likelihood  of  the  data 
under  the  continuous  inspection  Type  II  plan  for  the 
same  choice  of  r,  6 ,  and  .  The  efficiency  of  the 
interval-censored  Type  II  plan  is  computed  by 
computing  the  ratio, 

variance  under  continuous  inspection  plan 
mean  square  error  under  interval  -  censored  plan 
A  part  of  our  computational  effort  is  presented 
below. 


Sample  size  n=20 


r 

to 

e 

Efficiency 

10 

0.5 

1 

0.9878 

10 

0.6 

1 

1.0551 

10 

0.7 

1 

1.1103 

10 

0.8 

1 

1.1252 

10 

0.9 

1 

1.0882 

10 

1.0 

1 

1.0904 

15 

0.5 

1 

0.9537 

15 

0.6 

1 

0.9306 

15 

0.7 

1 

0.9502 

15 

0.8 

1 

0.9298 

15 

0.9 

1 

0.9116 

15 

1.0 

1 

0.8532 

4.  Conclusion: 

Even  for  a  moderate  sample  size  like 
n  =  20  and  for  a  moderate  value  of  r  =  10 ,  the 
interval-censored  Type  II  plan  is  as  good  as  the  Type 
II  plan,  if  not  better.  More  work  is  needed  to  make 
an  overall  recommendation. 

Acknowledgements:  The  authors  are  grateful  to 
Mr.  Jerry  Stockrahm  for  a  careful  reading  of  the 
paper. 


3.  Simulations:  The  mean  square  error  of  the 
maximum  likelihood  estimator  of  6  under  interval 
censored  Type  II  plan  is  intractable.  We  evaluated 
the  mean  square  error  of  the  maximum  likelihood 
estimator  empirically  by  mounting  Monte  Carlo 
studies.  The  inputs  arc: 

n  =  sample  size  (10,20,30,40); 
r  =  no.  of  failures  allowed  to  be  observed; 
0=0.1,1,10; 

^0  =  length  of  the  time  interval. 

For  each  choice  of  r,  0 ,  and  ,  5,000  samples 
each  of  size  n  were  drawn,  maximum  likelihood 
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L  Introduction;  The  Type  II  plan  is  used  on  many 
occasions  in  order  to  estimate  the  lifetime 
distribution  of  a  product  under  investigation.  Begin 
the  process  by  choosing  two  positive  integers  r<  n. 

Select  a  sample  of  n  units  of  the  product,  set  them  to 
work,  and  observe  the  units  continuously  until  r 
units  fail.  The  objective  is  to  estimate  the  lifetime 
distribution  of  the  product  using  the  data  on  r  failure 
times. 

We  want  to  offer  a  modification  of  this  plan  in 
response  to  a  past  consultation  problem.  This 
problem  arose  from  two  diverse  fields:  one  from 
engineering  and  the  other  from  ornithology.  In 
order  to  expedite  observation  of  the  r  units  failure 
times,  only  periodic  inspections  were  made.  In  such 
a  plan,  the  exact  failure  times  of  the  r  units  will  be 
unknown,  but  we  will  know  how  many  units  failed 
between  each  of  the  consecutive  inspection  times. 

Formally,  the  inspection  plan  can  be  described 
as  follows.  Choose  and  fix  a  number  /q  ^  ^  •  Select 
a  sample  of  n  units  and  set  them  to  work.  Inspect 
the  units  at  times  /q  ,2/o  ,  ■  •  •  until  r  units  fail.  Let 
M  denote  the  number  of  inspections  needed.  Let 
=  Number  of  units  failed  during  the  i^^  inspection 

interval  ((/ - 

/=1,2,3 .  The  data  consist  of 

M,  Xi,  Xj,  X}^. 

These  random  variables  satisfy  the  following 
conditions; 

X\  +  X2+...+Xj^_x  <r  -\, 

and 

X\  +  X2-^.  .+Xf^  >r. 

We  call  this  plan  Interval-Censored  Type  II  plan. 

The  basic  questions  we  were  asked  to  address  were: 

1.  Evaluate  the  loss  of  information  in  some 
meaningful  way  if  one  adopts  the  Interval-Censored 
Type  II  plan  over  the  traditional  Type  II  plan; 

2.  Provide  some  guidelines  as  to  the  choice  of  /q  . 


In  this  paper,  we  attempt  to  answer  these 
questions  to  the  best  of  our  ability.  This  is  ongoing 
work  and  we  want  to  report  what  we  have  achieved 
so  far.  Let  T  be  the  underlying  lifetime  variable 
associated  with  the  product.  We  assume  that  T  has 
an  exponential  distribution  with  unknown  parameter 
0>Q.  The  probability  density  function  of  T  is  given 
by 

/^(x)  =  x>0,  e>Q. 

We  obtain  the  likelihood  estimator  of  0  under  the 
Interval-Censored  Type  II  plan  and  compare  its 
variance  with  the  variance  of  the  unbiased  estimator 
of  0  built  on  the  maximum  likelihood  estimator  of 
0  under  the  Type  II  plan.  This  comparison  is  done 
using  extensive  simulation  studies.  Some  relevant 
distribution  theory  and  asymptotics  are  presented 
here. 

2.  Some  Distribution  Theory;  Let 
be  iid  copies  of  T  and  ^1)  <  ^2)  <•••<  T'in)  the 
corresponding  order  statistics.  Under  the  Type  II 
plan,  the  data  consist  of  7(j),  7(2),  •••,  T(^r)- 
maximum  likelihood  estimator  of  0  is  given  by 

Under  6,  2r6l  6 has  a  chi-squared  distribution 

A 

with  2r  degree  of  freedom.  The  estimator  0^^  is 
biased,  but  the  bias  is  correctable.  More  precisely, 

=  (''/('■  -  0)^for  all  ^  >  0. 

The  unbiased  estimator  ((/"  —  has 

variance  9^  lir  -  2)  provided  r  >2.  (Epstein 
and  Sobel  (1953).) 

Let  us  focus  on  the  Interval-Censored  Type 
II  plan.  We  need  to  determine  the  joint  distribution 
ofthedata  M,  Xy,  X2,  ■■■,  Xj^.  We  proceed  in 
two  stages.  First,  we  obtain  the  distribution  of  A/ 
and  then  the  conditional  distribution  of 
. X„|A/  =  m. 


Distribution  of  M:  VXq{M  =  1)  =  Pr^(^l>  0 
=  Pr^  (at  least  r  failures  in  the  interval  (0,  /q  ]) 
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(1965).)  For  m>2, 
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Z  Xi>r.  For  the  sake  of  simplicity,  let 
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We  now  derive  the  maximum  likelihood 
estimate  of  0  based  on  the  data:  number  of 
inspections  made  and  the  number  of  units  failed  in 
each  inspection  interval,  i.e., 

M  =  m,X^  =Xi,X2  =X2,...,X^  =x„. 

The  likelihood  L  of  the  data  is: 

L  =  Constant(l  - . . 

(n.  \ 
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The  log  likelihood  is  given  by: 

In  L  =  Constant  +  ^  Z  j  j 
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v,=i  / 

The  derivative  of  the  log  likelihood  is  set  equal  to 
zero  in  order  to  obtain  the  maximum  likelihood 
estimate.  The  following  equations  achieve  the 
objective. 
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The  maximum  likelihood  estimator  of  0 ,  upon 
replacing  the  data  by  the  corresponding  random 
variables,  is  given  by 
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The  next  objective  is  to  obtain  the  asymptotic 
variance  of  6 .  Rewrite  the  derivative  of  the  log 
likelihood  as 
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The  asymptotic  variance  of  9  is  given  by 
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The  formula  for  the  asymptotic  variance  simplifies  to 
evaluating  successfully  .  We  will 

evaluate  the  expectation  using  the  conditional 
expectation  argument.  Note  that 
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The  critical  step  is  the  evaluation  of  the  conditional 
expectation; 
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estimate  of  6  computed,  and  the  empirical  mean 
square  error  evaluated.  The  objective  is  to  compare 
the  mean  square  with  the  variance  of  the  unbiased 
estimator  of  0  based  on  the  likelihood  of  the  data 
under  the  continuous  inspection  Type  II  plan  for  the 
same  choice  of  r,  0 ,  and  .  The  efficiency  of  the 
interval-censored  Type  II  plan  is  computed  by 
computing  the  ratio, 

variance  under  continuous  inspection  plan 
mean  square  error  under  interval  -  censored  plan 
A  part  of  our  computational  effort  is  presented 
below. 


r 

Sample  size  n-10 

e 

Efficiency 

10 

0.5 

1 

0.9878 

10 

0.6 

1 

1.0551 

10 

0.7 

I 

1.1103 

10 

0.8 

1 

1.1252 

10 

0.9 

1 

1.0882 

10 

1.0 

1 

1.0904 

15 

0.5 

1 

0.9537 

15 

0.6 

1 

0.9306 

15 

0.7 

1 

0.9502 

15 

0.8 

1 

0.9298 

15 

0.9 

1 

0.9116 

15 

1.0 

1 

0.8532 

4.  Conclusion: 

Even  for  a  moderate  sample  size  like 
n-10  and  for  a  moderate  value  of  r  =  10 ,  the 
interval-censored  Type  II  plan  is  as  good  as  the  Type 
II  plan,  if  not  better.  More  work  is  needed  to  make 
an  overall  recommendation. 
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3,  Simulations:  The  mean  square  error  of  the 
maximum  likelihood  estimator  of  0  under  interval- 
censored  Type  II  plan  is  intractable.  We  evaluated 
the  mean  square  error  of  the  maximum  likelihood 
estimator  empirically  by  mounting  Monte  Carlo 
studies.  The  inputs  are: 

n  =  sample  size  (10,20,30,40); 

r  =  no.  of  failures  allowed  to  be  observed; 

(9=0.1,1,10; 

(q  =  length  of  the  time  interval. , 

For  each  choice  of  r,  6 ,  and  /q  ,  5,000  samples 
each  of  size  n  were  drawn,  maximum  likelihood 
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